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Abstract 

Consider a transportation problem with sets of sources and sinks. There are profits 
and prices on the edges. The goal is to maximize the profit while meeting the following 
constraints; the total flow going out of a source must not exceed its capacity and the total 
price of the incoming flow on a sink must not exceed its budget. This problem is closely 
related to the generalized flow problem. 

We propose an auction based primal dual approximation algorithm to solve the problem. 
The complexity is 0{e~^{n^ + n log m)m log J7) where n is the number of sources, m is the 
number of sinks, U is the ratio of the maximum profit/price to the minimum profit/price. 

We also show how to generalize the scheme to solve a more general version of the problem, 
where there are edge capacities and/or the profit function is concave and piecewise linear. 
The complexity of the algorithm depends on the number of linear segments, termed C, of 
the profit function. 

1 Introduction 

The Transportation Problem is a fundamental problem in Computer Science. It was initially 
formulated to represent the problem of transporting goods from warehouses to customers. The 
mincost version of the problem is also known as the Hitchcock problem after one of its early 
formulators |15] . It is well-known that while the Transportation problem can be seen as a special 
case of the mincost flow problem, the latter can also be transformed into the former |21j . 

The first primal-dual algorithm, called the Hungarian method [T7|, was also proposed for 
the Assignment Problem which happens to be a special case of the Transportation Problem. 
More recently, the Transportation Problem was used by Bertsekas [3| and Goldberg and Tarjan 
(llj as a natural framework for applying auction based algorithms. 

The typical transportation model assumes that there is no loss of goods as the goods are 
transported from the sources to the sinks. In real life, there is often a loss involved. For example, 
in case of electricity there is loss due to resistance and in case of oil, due to evaporation. In case 
of a delicate merchandise, there might breakage. Sometimes, there may even be a gain (e.g. 
transfer of money: currency conversion). The loss or gain can be modeled as a price function 
on the transport from the source to the sink. 

As an example, consider a set of depots and a set of retail stations. The depots have a 
certain supply of goods, while the retail stations have a limit on their intake capacity. For every 
unit transported between a pair of depot and station, a certain amount of profit is made which 
depends on selling price, cost of transportation etc. In addition, consider the situation where 
goods transported suffer a loss proportional to the number of units transported. This factor is 
unique to each pair as it depends on the mode of transfer and distance etc. In order to measure 
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the net incoming amount at the stations, we have to account for this loss. The intake capacity 
hmit is therefore, on a weighted sum rather than a simple sum. 

In a more recent context, keyword based advertising is widely applied by internet based 
service companies like Google and Yahoo. Consider the following model for assigning adver- 
tisement space to the bidders: each bidder, say j, modeled by a sink, specifies its total budget 
bj and the price pij it is willing to pay for the ith keyword, specified by source i. An estimate 
of the number of times a keyword would be invoked in web searches provides a capacity Oj for 
the ith keyword. Assignment of a keyword i to a bidder j provides a profit Cij. The goal is to 
maximize the revenue of the service provider without exceeding the bidders budget. A version 
of the problem restricted to an integral setting may be found in [2] . The online version of the 
problem has been considered in |19| . 

In order to model the above and similar problems, we propose a generalization of the trans- 
portation problem: the Budgeted Transportation Problem (BTP) which is closely related to 
Generalized Flow Problem. We then present an auction based approximation algorithm. We 
also show how to extend the technique to further generalizations of the problem. 

Like the conventional Transportation Problem, the problem is modeled on a weighted bi- 
partite graph. In our version of the problem we consider maximizing the profit in transport- 
ing/assigning goods from the sources to the sinks, subject to the following conditions: Each 
source, i has a supply of goods, bounded by a capacity function Oj, and each sink j has a bound 
on the incoming set of goods, specified by a budget bj. The difference from the conventional 
transportation problem is that the budget is a bound on the weighted sum of the incomming 
goods as follows: each source-sink pair has a price function pij and the bound bj is on 

the sum YliPijfij where fij is the amount of goods transported/assigned from source i to sink 
j. The goal is to maximize the sum of profits of transporting/assigning goods from sources to 
sinks, given that the profit obtained by transporting/assigning fij units is Cijfij. The problem 
can be formally stated as follows: 

Given a bi-partite graph G{S,T), capacity a : 5 — >■ budget b : T ^ Z+, price p : 
5 X T — )• Z^, profit c : 5" X T — )■ Z"*" the budgeted transportation problem is to find a fiow 
function f : S x T ^ so as to maximize YlieS jeT'^ijfij subject to capacity constraints, 
Ylj^T fij — ^^'^ budget constraints YliesPijfij — ^ji^i ^ T.We let |5| = n and |T| = m. 

Note that, for simplicity, the price function, like the other functions, is assumed to be integral 
in the paper. This could be replaced by a price function to the set of rationals. 

We describe an approximation scheme for the budgeted transportation problem. Our algo- 
rithm is a primal-dual auction algorithm. Auction algorithms have been utilized before in the 
solution of combinatorial problems |3l [11] including the classical transportation problem [5j, 
and more recently for finding market equilibrium [9l [10] . Bertsekas et al ^ apply the auction 
mechanism to solve the generalized fiow problem but the complexity of that algorithm is not 
polynomial. 

While the auction technique has been used before, the application in the current context is 
different. The auction mechanism is used to realize a set of tight edges, i.e. edges which satisfy 
complementary slackness conditions. Paths and cycle are found in a subgraph in these set of 
tight edges. The fiow path/cycle could either be fiow generating or create a fiow reduction. The 
primal-dual nature of our approach allows us to push fiow along these paths or cycle without 
worrying about the profit of the path/cycle. The cycle can be eliminated in linear time. Our 
algorithm may be interpreted to have achieved a version of cycle-canceling [2S] without having 
to compute the associated profit. 

We achieve a complexity of 0(e~"'^(n^ -|- nlogm)mlog t/logm), where n is the number of 
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maxij(^) 

sources, m is the number of sinks and U is . The problem can also be modeled as 

emmij(-J-) 

a generalized flow problem. This can be done by adding a super source connected to all the 
sources and a supersink connected to all the sinks and minimizing the negative of the profit 
function. By using the technique of scaling as in [18J, the dependancy on \ogU can be replaced 
by log(nm-/e). We will thus choose to ignore terms involving logf/, and in fact logn, when 
comparing our results with other results. 

There are numerous results for maximum/mincost generalized flow and related problems 
[21 E EHl Ea EOl [la [13]. Some of the currently best known FPTAS for the generalized flow 
problems are of complexity 0{\oge~^E{E + Vlogl))) [3 [2l [23] for maximum generalized 
flow and 0{e~'^E'^VJ) [7] and 0{loge^^ E'^V'^) [25j for minimum cost generalized flow. Here, 
/ = logM, J = log log M + loge^^ and M is the largest integer in cost representation. For 
the special case, when there are no flow generating cycles, the complexities are 0(e~^i?^) and 
0{e~'^E^J) for maximum and minimum cost generalized flows respectively [3, where V is the 
number of nodes and E is the number of edges in the given graph. 

A straight-forward application of the best known generalized flow approximation algorithm 
to this problem would require 0(log e~^{mn)'^{m+n)'^) time using the minimum cost generalized 
flow algorithm in [25] or 0{e~'^{mn)'^{m + n)J) time using the algorithm in [7J. Our algorithm 
has a better complexity than the first when e > and is better than the second for all 

ranges of e. 

If we use the packing algorithm from [8j than one can solve this problem in 0(e^^(n + 
m)nm) since we have m + n rows and mn columns in the primal constraint matrix. Other 
known combinatorial algorithms for packing |22[ are also all dependent on . After the 
completion of this work we became aware of another result [IB] which provides a randomized 
(1 + e) approximate solution to the packing problem in time 0{N + (r + c) log(A^) / e^) where N 
is the number of non-zero entries in the specification of the packing program and r + c is the 
number of rows and columns in the matrix specifying the packing program. As applied to the 
budgeted transportation problem this algorithm would require 0{nm{lognm)/e^) and would 
be faster than our approach when e > 1/n. 

Comparing with algorithms that have a dependence on 1/e, one recent algorithm for packing 
problems is by Bienstock and Iyengar [6], which reduces the packing problem to solve a sequence 
of quadratic programming which in turn are approximated by linear programs. Consequently, 
the complexity will have high polynomial dependance on n and m. We achieve a combinatorial 
algorithm with complexity dependant on 0(l/e) and a low polynomial in n and m. 

We extend our algorithm to solve a generalization of BTP called the Budgeted Transship- 
ment Problem (BTS). It is similar to BTP except for added edge capacity constraints. We 
achieve the same approximation with no additional time complexity, i.e 0{e~^{n'^+n log m)m log C/ log m). 
This is of improved efficiency, since, again, to compare this with existing results: the result in 
|25j would imply a complexity of 0(log e~^(?Tin)^(m -|- n)^), [2| would imply a complexity of 
0(e~^(mn)^(m -|- n)J) while [8] would imply a complexity of 0(e~^(n + m -\- mn)nm) as the 
number of constraints is now m + n + mn. The results from [16] would, again, be of complexity 
0(nm(log nm) / e^). 

In section [XT] we present a basic auction based algorithm and show its correctness. In section 
12.21 we describe a modified algoritm and discuss its convergence. Since the algorithm is similar 
to the algorithm for the generalized problem, BTS, we demonstrate an outline of the algorithm. 
Section 2] extends the result to BTS. In section 14.11 we present a detailed algorithm for BTS 
(which also applies to BTP) and in subsequent sections the proofs of correctness and complexity 
of the algorithm. The piece-wise linear case is discussed in section [5l We conclude in section [6l 
We also show a relationship between BTP and minimum cost generalized flow problems in the 
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appendix (section [Tj) . 

2 Auction Based Algorithm 
2.1 The Basic Method 

In this section, we describe a primal-dual auction framework to approximate the maximum 
profit budgeted transportation problem (BTP) to within a factor of (1 — e). The algorithm is 
a locally greedy algorithm where the sources compete against each other for sinks by bidding 
with lower effective profit, a quantity defined using dual variables. This algorithm will serve to 
illustrate the methodology. In order to prove complexity bounds we will, in the next section, 
refine the algorithm further. 

We first model the problem using an LP and also define its dual. 



maximize (^ijfij 

i€S,jeT 

subject to : 

fij <ai yi e S 

^Ptjfij<bj ViGT 

In the remainder of the paper, the variables i and j will refer to a source i in S" and sink j in 
T, respectively. 

The dual of the above program is 

minimize ajOj + /3jbj 

subject to : 

> Cij - pij/3j yij 
ai,f3j > 

The following variables based on primal and dual variables are used in the algorithm. 

• Oj : the dual variable associated with each source. 

• /3j : the dual variable corresponding to the budget constraint on the sink. An algorithmic 
interpretation of /3j is that it represents the value of sink j. It increases as different sources 
compete in order to ship flow to it. At every instance of the algorithm, flow may be pushed 
onto the sink j from a source at a given value of 13 j. We term this an assignment of the 
source at valuation /3j. Further, at any instance of the algorithm sources are assigned to 
sinks at two prices (values), some at value P while others are assigned at a companion 
valuation /Sj where /3j = /3^(1 + e). 

• Uij: the valuation of sink j when assigned to source i to ship the flow from i to j. It 
equals /3j when source i ships flow to sink j. Subsequently, as /3j rises, it becomes equal 



4 



• Si : the supply (or surplus) remaining at the source i. Sj is evaluated as Si = ai — fij 
and is implicitly updated in the algorithm. 

• dj : the amount remaining from the budget at sink j. dj is evaluated as dj = bj — Y2iPijfij- 

Because the input is asssumed to be integral, we note that all the variables are rationals, 
with a bound on the absolute value of the denominators. 

The initial and primal dual feasible solutions are obtained by the initialization procedure [2TT] 
Initialize. The algorithm maintains a pair of feasible dual and primal solutions. It works to 
modify them until the following complementary slackness conditions, which define optimality, 
are exactly or approximately satisfied. 

CSl : fijicij - pijPj - Oj) = Mij (1) 
CS2: ai{ai-Y.jfi3) = ^ (2) 

CSS: /3j(6i-EiPii/^,) = Vi (3) 

In each iteration, the algorithm 12.21 Auction picks a source, i, with positive surplus, Sj = 
Qi — fij, and Oj > 0. It then picks a sink j that provides it the largest effective profit, where 
the profit is specified by the quantity, Cij — pij/3j. If the sink is not saturated flow is pushed 
from source i to sink j. If sink j is saturated, flow from a source i' , that is shipping flow to 
sink j at a lower of value f3j, is replaced by flow from i, i.e. source i wins over the shipment to 
sink j from source i' by bidding a higher value (/3j) for j. After the change of flow assignment, 
it is checked (|2.3I UpdateBeta( /3j ) ) that there exists at least one source shipping flow at 
the lower value. If not, then /3j is raised, indicating that shipping flow to j requires paying a 
higher value (and lower effective profit). At this stage, the process is repeated considering sinks 
in order of effective profit until 
(i) either all the surplus at node i is pushed out or 
(n) ai = 

Algorithm 2.1 Initialize 

1: Oi maxj Cij 

2: f^j = 
3: /3j ^ 
4: ^j ^ 



To prove the correctness, we show that primal and dual equations are satisfied and com- 
plentary slackness conditions hold. All line numbers in the discussion below correspond to the 
description of Algorithm 2.2 Auction, except when stated. 

Lemma 2.1 The dual solution maintained by the algorithm is always feasible. 

Proof: We observe that the value of Pj never decreases as the algorithm progresses. Therefore, 
oti < Cij — Pij(3j can only become true when ai is modified. However, when Oi is modified, it is 
set to a value > maxjjcjj — Pij(3j} in line [231 At that point, Oi > Cij — Pij(3j for all j. □ 
For use in the next lemma, we define the total price of the incoming flow on sink j to be 

J2iPijfij- 

Lemma 2.2 The primal solution maintained by the algorithm does not exceed the budget con- 
straints. 
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Algorithm 2.2 Auction 



while there exists a source i such that > and fij < ai do 
Pick a sink j such that Cij — Pijf3j is maximized 
if the sink is saturated, i.e. ^iPijfij = bj then 

Find a source i', currently assigned at the lower value level i.e. yi/j = Pj 
ii i' ^ i then 

//determine the amount that can be replaced 
X ^ min(si, UjPi'j/pij) 
/ / replace the flow 

fij fij ~^ 

fi'j ^ fi'j -^Pij/Pi'j 

Vij ^ Pj 
else 

/ /i is already shipping to j at the lowest level 
//Raise the assigned value without changing any flow 

Vij ^ Pj 
end if 

Update Pj and Pj 
else 

//The sink is unsaturated, push maximum flow possible 
//under supply and demand constraint 
X ^ mm{si,dj/pij) 

fij fij "I" ^ 

end if 

ai maxj{cij - Pijl3j,0) 
if = then 
Vij ^ Pj for all j 

//This ensures there is no further rise in without reducing fij to zero 
end if 
end while 



Algorithm 2.3 UpdateBeta {Pj ) 

1: if j3j = then 

2: Pj emin, Cij/pij 

3: else if Vi : fij > 0, yij = f3j then 

4: /3j ^ /3, 

5: Pj^Pjil + e) 

6: else 

7: //no update required 
8: end if 
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Proof: The flow on any edge is increased in lines [9] and [22l In line [9] the total price of the 
incoming flow on sink j increases by xpij. It is then reduced by piijX^^ in line llOl and therefore, 

the budget is met. 

In lineEIl the total price increases by at most dj = bj — J2iPijfij thus the budget is 
not exceeded. □ 

Lemma 2.3 The primal solution maintained by the algorithm does not violate supply con- 
straints. 

Proof: During the algorithm, the increase x in flow, fij, on an edge ij is always limited by the 
available supply (line 7 and 21) and therefore, the supply constraint is clearly satisfied. □ 

Combining the above lemmas |2 . 2 1 and 1 2 . 3 1 we conclude: 

Lemma 2.4 The primal solution maintained by the algorithm is always feasible. 

Lemma 2.5 During the course of the algorithm, if sink j is not saturated then (3j = 0. Further, 
Vj, j3j is non- decreasing during the course of the algorithm. 

Proof: The variable Pj for all j is initialized to zero. It is subsequently changed by the pro- 
cedure UpdateBeta, which in turn is called only if the sink j is saturated. That is, line [T7| is 
executed only when the "if condition in line [3] : YliPijfij — is satisfied. Moreover, Vj, /3j 
only increases in value. □ 

The following lemma will be useful for bounding the complexity of the various algorithms: 

Lemma 2.6 Suppose fij, the flow on edge ij, that was pushed at a valuation of sink j, yij = 
Pj = + e), is reduced to zero, i.e. pushed back to zero. Then at the next occurrence of the 

event that corresponds to fij being pushed back on edge ij, the valuation of sink j, yij > 

Proof: This follows from the fact the /3j is monotonically non-decreasing. Flow is pushed back 
on an edge ij only when there is flow on the edge, fij that had been pushed at the companion 
valuation of Pj/{1 + e), where (3j is the current valuation. Let flow be pushed back from sink 
j to source i on edge ij at some time t in the algorithm such that the flow is reduced to zero. 
Then when flow is pushed from i to j, at a subsequent time t' > t, the valuation yij > fij (line 
11), since the valuation of j does not decrease. Consequently, if flow is pushed back again from 
j to i, the valuation has increased. □ 

Lemma 2.7 During the execution of the algorithm, the Complementary Slackness condition, 
CSl is approximately satisfied for each edge ij. That is, 

^ij ■ fij > 0, < Cij - pijPj + ecij 

Proof: When flow is pushed on the edge ij, fij > and Oi = Cij — PijPj or is (line 2 and 
24 in Algorithm 2.2 Auction) and the condition is true. Note that at this event, yij = /3j. 
Subsequently, Oi or Pj may change. We consider the changes to and Pj and to the flow, 
fij, between the event in the algorithm when positive flow is pushed on the edge ij and the 
event when this pushed flow reduces to zero, being pushed back on edge ij. In between the two 
events, since we assume that fij > 0, yij is either /3j or /3j/(l + e). We consider cases depending 
on whether Oj > or not. 
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1. Oj > 0. Suppose ai does not change but f3j has increased by a factor of (1 + e). In this 
case Uij has a value equal to Pj, otherwise fij would be zero. Since the value of /3j changes 
by a factor of (1 + e), < Cij — Pijf3j + ecij. Note that the value of /3j can not increase 
further unless fij reduces to zero. 

If Oj has also changed, it only decreases, since /3j is increasing. Therefore, the inequality 
still holds. 

2. Oj = 0. No future change in the value of Oj can occur. Prior to the step at which is set 
to zero, ai = Cij — Pij/3j > for some j. This implies that Pij(3j < Cij. An increase in the 
value of I3j thus still results in Cij — Pij(3j + ecij > 0. Thus the inequality holds when ai is 
set to zero. Further, when ai becomes zero, the algorithm sets the value of yij to be /3^- so 
as to ensure that Pj is not increased without the event that flow is pushed back on ij and 
fij is reduced to zero. And, subsequently, source i is never considered for pushing flow 
fij since line 1 of Algorithm 2.2 Auction only considers vertices with ai > 0. Thus if 
flow fij > then (3j has not changed. And if f3j has changed, the flow is set to zero. The 
inequality thus holds. 

□ 

Furthermore, 
Lemma 2.8 Algorithm 2.2 Auction terminates. 

Proof: At each iteration of the algorithm, flow is pushed on an edge and Vj, f3j either remains 
the same or the value of (3 increases for some j (Lemma [23]). 

When flow is pushed on an edge zj, either (i) Si goes to zero (ii) the budget at a sink dj is 
met or (iii) the flow on a back edge ji' is reduced to zero. Note that, by lemma [231 in case (iii) 
flow cannot be pushed (in the forward direction) on the edge i' j without increasing the value 
of yij. 

The first event, case (i), results in flow being pushed from i to j such that yij = f3j and 
will happen only once per source-sink pair until an event where flow is pushed back on an edge 
j'i to add to surplus at i occurs. Thus this event occurs whenever new surplus is generated, 
which is bounded by the number of times flow is pushed back on an edge without case (iii) 
ocurring. The number of times flow is pushed back on an edge, without reducing the flow to 
zero (case (iii)), is bounded since the amount of change of flow at every instance is a rational 
with bounded denominator. 

The second event, case (ii), happens at most 0{m) times, since after that event the sink 
remains saturated during the subsequent operations of the algorithm. The third event, case 
(iii), occurs at most n times before all edges leading to a sink have flow reduced to zero, and 
consequently the value of Pj must rise. This bounds the total number of operations in the case 
that j3j remains the same. Further, the number of increases of f3j is bounded since f3j increases 
by a factor of l + e and its value is bounded above by maxj{cjj}. Thus the algorithm terminates 
in a finite number of steps. □ 

In order to show good convergence of the algorithm we introduce further modifications. 
However, when the algorithm terminates then the following is true: 

Lemma 2.9 At termination, the difference between the value of primal solution ^ijCijfij and 
the value of the dual solution is at most e^^ - (^ijfij 
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Proof: The value of the dual solution is 

= ^mai + ^bj^j (4) 

i j 

i j j i ij ij 

When we subtract the value of the primal solution Cijfij, from the above, the difference 

is 

'^iat - fij)ai + '^{bj - Ypijfij)Pj - X] fij{cij -ai- pijPj) (6) 

i j j i ij 

Thus, the total absolute difference is at most Ai + A2 + A3 where 

^1 = I X] -^^J'^^^J ~ ~ Pijl^j)\ 
ij 

A2 = l^Oiiai - 

i j 

^3 = \Y(3,ibj-Yp^jfij)\ 
j i 

From lemma [221 we have 

Vij : fij > 0, Qi < Cij - pijf3j + ecij 

Therefore, 

Ai < e^^Cijfij 

ij 

Prom the termination condition of the algorithm, we know that for any source i such that 
«i - J2j fij > 0' we have ai = 0. 
Therefore 

A2 = 

For all unsaturated sinks f3j = 0. Therefore, 

A3 = 

Combining all of the above we get the lemma. □ 



2.2 Modified Algorithm 

In order to show the required complexity bound, and make its analysis simpler, we make a few 
modifications to the above algorithm. We address the details of the data structures required 
and present the complexity analysis. We also show that the modifications do not affect the 
correctness of algorithm 12.21 Auction as proved in section [21 

We introduce the concept of a preferred edge, back edge and derived graph as explained 
below. 
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Preferred Edge: For a given set of values for and Pj we designate a preferred edge Pri 
for each source i. This edge maximizes Cij —pijPj, i.e. Pr^ = ik where k = argmaXj(cjj —pij/3j). 
In the case that multiple values of j satisfy the condition, only one is picked. The source i will 
be required to push flow along this edge until the edge is no longer is preferred. 

Back-edge: All edges ij such that fij > and yij < fij. Flow is pushed back on these edges. 
The set of all such edges is termed B. 

Derived Graph: We define the capacitated derived graph H with respect to an intermediate 
solution maintained by the modified algorithm. This graph consists of all the sources and sinks 
in the given bipartite graph and directed edges which are either the preferred edges or the back 
edges. All the preferred edges are oriented in the forward direction, whereas all the back-edges 
are oriented in the reverse direction, indicating the directions along which flow will be pushed. 
The capacity of a preferred edges is infinite whereas the capacity of a back-edge is the amount 
of flow carried on that particular edge. Formally, H = (S, T, E, c) where 

E = {ij\i G 5, j G T, ij G Pr,} U {ji,j eT,ieS, ij G B} 

and the capacity function c : E Z"^ is as follows: 

c{ij) = oo iff ij G Pri 

c{ji) = fij iff ij G B 

Before starting with the main algorithm, the following preprocessing step eases the descrip- 
tion: 

Pre-processing step: Consider a preferred edge ij G Pvi. If edge ij is also a back edge (thus 
creating a 2-cycle {i j ^ i) in H) and is not the only back edge from j then the valuation 
at which flow is shipped from source i to sink j, i.e. yij = fij /{I + e), is increased to fij. Note 
that this modification is not done when the sink has only one back edge ij since the existence 
of at least one valuation equal to /3j/{l + e) does not allow a raise in the price of sink j. This 
ensures: 

1. If there is a 2-cycle on sink j, sink j has only one back edge. 

2. /3j does not change. 

In the derived graph H, if ij G Pri before the modification then ij G Pr, after the modification 
also, and source i will choose to push fiow along its preferred edge. 

The algorithm in section 12.11 pushed flow from a source to a sink causing flow to be pushed 
back to some other source and consequent independant processing of that other source was 
performed. The modified algorithm in this section, starting from a source with a surplus and 
positive ai, finds a path along the edges of the derived graph H and pushes the fiow according 
to the type of path discovered, thus ensuring that the surplus fiow is processed in one push 
of fiow along a fiow path. The path found by the algorithm ends at either (i) an unsaturated 
sink,(ii) a source with = or (iii) a cycle. 

We detail the various cases that arise. A source i is picked such that Oj > and fij < ai. 
And a path, P, is discovered by, say a depth first search, along the edges of the derived graph 
H. There are different types of paths/cycles that may arise. 



10 





Figure 2: An example where path ends in a sink with Pi = 
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Figure 3: Example where the path ends in a cycle 



1. Type-1 paths: We have two cases, one where the path finds a source where = 
and the second where the path finds an unsaturated sink jk where =0. In the first 
case, let P = io, jo, . . . and in the second case let P = io,jo,---,jk- In both cases 
flow is pushed as follows: in the first case flow is pushed into source i/^ (i.e. returned 
to the source) and in the second case flow is pushed into a sink (jk), without violating 
feasibility or complementary slackness conditions. However, the flow may be limited by 
the capacity of the edges on the path, P. Wc define the capacity of the path, c(P) as 
c(P) = min{c(e)|e G i? fl P} The processing of the flow is as follows: 

(a) The flow is limited by the surplus, at the starting source zq. Push all the surplus 
along the path. The surplus at the source disappears. This is true for both cases. 

(b) The amount of flow is limited by the capacity of an unsaturated sink on path P 
(second case), i.e. c(P) > dj^,. In this case, we push enough flow to saturate the sink 
and consequently, raise the value of Pj. 

(c) The amount of flow is limited by the capacity of the path c{P). In this case, we 
push flow from the source along the path. Let / be the flow starting at the source. 
Let e = jiii+i be the first edge such that / > c(e) = c(j;,i;+i). The amount of flow 
that can be pushed along this edge is limited by c(e). Pushing flow along this edge 
reduces the capacity of the edge e to 0. The surplus / — c(e) is added to the surplus 
of the previous source i; which has iiji as the preferred edge. The flow pushed back 
through jiii+i is accumulated at ii+i as a surplus. The surplus at i^+i is pushed 
further along the path recursively. In the end one or more back edges would have 
zero flow and surplus may be left at each source ocurring just prior to the sink where 
the back edge with zero flow is incident. 

2. Type-II paths: The path ends at a 2-cycle A 2-cycle implies ji is the only 
back-edge on sink j, as soon as the flow is again pushed along edge ij, yij changes (line 
15, Algorithm 2.2 Auction) and the only back edge disappears. This results in an 
increase in the value of fij. 
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3. Type-III paths: A cycle is discovered since one of the sources is encountered again. Let 
the path discovered be denoted by iq, jo.-.ik---jqik which is decomposed into a simple cycle 
C = ikjk ■ ■ ■ jqik aiid a simple path Pi = io, jo ■ ■ ■ ik- The source ik being the entry point, 
i.e., the first source that appears in the cycle. First the flow is pushed along the path 
io...ik taking into account the capacities as in Type-I paths. Any additional surplus that 
appears at source is handled as follows. 

Consider an imaginary procedure that simulates sending flow around the cycle, C, multiple 
times, while conserving the flow at each of the intermediate sources in the cycle. 

We define the transfer ratio pi for a back-edge w.r.t. C to be PiiJi/Pii+iji where iiji 

is the preferred edge for source i and jiii+i is a back edge. 

Further define the cumulative transfer ratio pf for the back-edge ji,ii^i with respect to a 
cycle C to be the amount that goes through the edge ji,ii+i if one unit of flow is pushed 
along the cycle starting at source ij^. This is equal to Y[^z=k Pz- 

Let the cumulative transfer ratio of the cycle be Pq = H^La; Pi ■ This is the amount of flow 
that reaches the first source, ik, back after being pushed all around the cycle starting with 
a unit of flow at ik- 

Let us call the push of flow once around the cycle as one revolution. If we were to send 
the surplus at source ik, Si^ units of flow, repeatedly in this cycle for say, r times, while 
conserving the flow at each of the intermediate sources, the total flow shipped through 
the back edge jiii+i is given by the geometric series F{1, r) = Yl^z=Q ^ikPiiPoT- Since this 
back edge is capacitated, the maximum number of times flow can be pushed around the 
cycle before saturating this edge to capacity, which we term as the limiting number of 
revolutions, is given by Ri = argmaXj.(F(Z, r) < fii_^^ji)- This number can be computed 
in Oil) time. It ignores that other edges may also limit the flow. We thus compute 
Rmin = miii^gsnc Rl to compute the limiting number of revolutions for the cycle, i.e. the 
smallest value of the limiting number of revolutions, computed from amongst the back 
edges in the cycle, C. 

The flow on the edges of the cycle is modified as: 

fidi ^ fiui + ^hPf-iiPeV 

j=0 

for forward edges iiji and 

fk+ih ^ fk+iJi ~ X] ^iiPliPQy 

3=0 

corresponding to back edges. 

We consider cases based on the value of p©: 

(a) Pq < 1. In this case we have a decreasing gain cycle, C, i.e. the value of the 
flow decreases as it is sent around the cycle. We compute the limiting number of 
revolutions for each back edge in the cycle, C. Note that this number may be infinity 
for some edges. However, this can be determined by computing the value to which 
the geometric series, F{1,), converges and does not require simulating an infinite 
sequence of pushes. 
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After Rmin rounds, the remaining surplus at is reduced to Si^pQ"^'". Note that if 
Rmin is infinity then the surplus at the source, Sj^, decreases to zero. The surplus at 
every other source in the cycle remains unchanged. 
Consider the cases depending on the value of the surplus at i^: 

(i) the surplus remaining at if^ becomes zero, or, 

(ii) the flow on one of the back edges reduces to zero. This case occurs when the 
number of revolutions is limited due to one of the back-edges. Since Rmin integral 
revolutions need not achieve this, one more sequence of flow pushes around the cycle 
may be required. 

(b) Po > 1. In this case the flow increases or remains the same as it goes along the cycle. 
There will be thus no decrease of surplus at the starting source vertex k. However, 
for each of the back-edges Ri = argmax^(F(/, r) < fii^^j,) is a finite number since 
each push of flow around the cycle reduces the capacity of the back-edges. The 
smallest of these values is chosen and the maximum number of rounds of flow to be 
pushed around the cycle is computed. Again, since Rmin integral revolutions need 
not achieve this, it may require one more sequence of flow pushes around the cycle 
to reduce the flow on a back-edge to zero. 

In the cases above where the surplus is not reduced to zero, at least one edge is saturated 
by the flow pushed around the cycle. Any surplus remaining at other nodes is pushed 
around the cycle to accumulate at the source(s) just before the saturated edge(s). 

The overall structure of the algorithm is similar to that in section [21 After pushing the flow 
through the path or cycle discovered, values of /3 are updated if required, preferred and back 
edges computed as necessary and the algorithm repeats until all surplus is removed or until 
complementary slackness conditions are satisfled. This algorithm is termed as the Modified 
Auction algorithm. 

Data Structures: In order to determine maxj Cij —pij/3j for each i we maintain a heap. This 
heap can be updated when /3j changes in O(logm) time. The UpdateBeta procedure, as in the 
previous algorithm, utilizes a set data structure to determine the need for update of (3j in 0(1) 
time. 

We next discuss how this approach improves the complexity. However, we can actually solve 
a more generalized problem, i.e. where each edge {u,v),u G S,v G T has a capacity n^. So we 
defer a formal description and proofs of the above approach to section [3l 

2.3 A discussion on the correctness and complexity of Modified Auction 

We discuss the correctness and time complexity of the modified algorithm. Formal proofs are 
provided when we consider the generalized model including edge capacities. 

Correctness: In order to justify that the modified algorithm terminates with a (1 — e) ap- 
proximate solution, first, we observe that while the algorithm 12.21 Auction does not require any 
particular order in which the surplus sources are picked, the sequence of choices of the sources 
from where the surplus is reduced, as imposed by the modifications, still ensures primal and 
dual constraints. 

In the case that the algorithm encounters cycles, the modified algorithm essentially simulates 
multiple steps of the algorithm 12.21 Auction. The end result is exactly the same if each step 
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was performed individually and repeatedly. The amount of flow shipped out of a source at 
each step, however, may be smaller than what Algorithm 12.21 Auction would have shipped out. 
However, the primal and dual constraints are still ensured. 

The modified algorithm makes the same changes to the variables ai, 13 j or fij as described 
in the previous section. The only modification is in the change to the variables, yij in the 
preprocessing procedure. This change does not affect the claim that f3j < yij{l + e) for all i in 
the proof of lemma [2771 as yij is set to (3j itself. The pre-processing procedure leaves at least one 
back-edge intact; this ensures that no change in (3j takes place as a result of the preprocessing. 
All the arguments involving f3j are, therefore, unaffected. Thus primal and dual feasibility is 
maintained. Further the complementary slackness conditions are satisfied. 

Complexity: For the complexity we consider the sequence of steps between two successive 
raises of the value of /?, termed as a phase, in the algorithm. 

We use the fact that for all j, the total number of times the value f3j rises is bounded. We 
show that each operation during the course of the algorithm can be charged to a rise in the 
value of /3j for some j. In fact, V/3j, every rise in the value of /3j is required to be charged 
at most 3n^ + n times during the algorithm. We account for the work required for each push 
of flow along a path or cycle via a charging argument. There are at most n pushes of flow 
along edges of a path before either: (i) a 2-cycle is encountered, (ii) a source with a = is 
encountered, (iii) an unsaturated sink is encountered, or (iv) a cycle is encountered. Note 
that in the last case all the calculations for determining the limiting number of revolutions can 
be done in 0(n) time, since the cumulative transfer ratio can be determined by traversing the 
cycle once. As noted above, in all the cases that arise when flow is pushed around the cycle, 
either the surplus that has been pushed into the cycle reduces to or the flow on one of the 
back edges is reduced to zero. The 0{n) amount of work can be charged to the disappearance 
of surplus at a source or reduction of flow to zero on a back-edge ji, and hence an increase in f3j, 
subsequently, when all back edges incident to the sink j have the corresponding flow through 
them reduced to zero. The other cases have a similar analysis. The detailed complexity of the 
algorithm will be described in the next section. 

3 Edge Capacities 

In this section we further generalize the Budgeted Transportation Problem. In this version each 
edge has an upper bound on the flow going through. The capacity of each edge is represented 
by Uij G Z"*". We call this problem the Budgeted Transshipment (BTS) problem. It is easy to 
see that the budgeted transportation is a special case where edge capacities are infinite. 
Consider the Linear Program for BTS 



maximize 




(7) 
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The dual of the above program is 

minimize ^ ajOj + ^ Pjbj + ^ Uijlij (12) 

i j ij 

subject to: 

ai > Cij - pijPj - -jij yij (13) 

ai,^,-,7y>0 yij (14) 

And the complementary slackness conditions are 

fij{cij -ai-pijf3j - -/ij) = Vij (15) 

<^i{(^i -^fij) = Vi (16) 
j 

^jibj-YlPi^fi^) = (1^) 

i 

lij{uij - fij) = Vij (18) 

4 Approximation Algorithm for BTS 

We now present a detailed approximation algorithm for the BTS problem. 
We first (re) define some terms 

Edge dual variable jij-. The dual variable associated with each edge. If ij is saturated, jij 
is set as, = minj/{(cy —PijPj) — (%' — Pij'Pj')} over all choices of edges if. It is set to if 
edge ij is unsaturated. 

Preferred Edge: For each source i an edge ij, termed Pri, such that ij is unsaturated and 

Cij — Pijf3j is maximized. Note that in the case when more than one edge meet the criteria, one 
of them is chosen as the preferred edge. The preferred edges are directed forward (from i to j) 
in the derived graph. Therefore, every forward edge is unsaturated. 
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Back Edge: An edge ij such that fij > and yij = j^. 

As in the previous algorithm, a derived graph H is used by the algorithm. The derived 
graph has vertices as the original vertex set and edges as back edges and preferred edges. The 
edges are directed with the preferred edge from source i, Pri = ij, being directed from source i 
to a sink j, with capacity given by Uij — fij, and a back edge directed from a sink j to a source 
i' with capacity /j/j. 

The algorithm differs from the one in section [2] in the use of the variables 'jij . To ensure 
complementary slackness, this variable becomes non-zero only if edge ij is saturated. This 
variable is not explicitly maintained but its value can be deduced from and Pj. Its value 
indicates the extra cost incurred in using another edge from source i. To ensure complementary 
slackness conditions, the value of the variable jij implicitly reduces when f3j rises and when the 
reduced value becomes zero, the edge ij can be designated as a back-edge such that flow can 
be pushed back on the edge ij, leaving it unsaturated. 

The following additions/modifications are made to the algorithm: 

• During the initialization, 7jj is set to zero. 

• f3j update: This update occurs when yij = f3j, for all unsaturated edges with positive flow. 
See procedure 14.41 update. 

• Preprocessing : Before the removal of 2-cycles, for each saturated edge ij, if jij > 0, ij 
cannot be a back edge. In case 'jij < 0, it becomes a back edge. (See procedure 14.41 
update). 

• Path/Cycle discovery and pushing of flow: In the BTS problem, the forward edges have 
a capacity as well. While pushing the flow in a path or in a cycle, the bottleneck edges 
could be forward edges in addition to back-edges. (See 14.71 findpath. 14.51 pushFlowPath 
and 14.61 pushFlow^CycIe) . 

The algorithm with the above modification is described as Algorithm 14.11 Modified Auc- 
tion, along with the associated procedures. 

The algorithm initializes the primal and dual solution using procedure 14.21 Initialize. Pre- 
processing is done in procedure 14.31 Preprocess to choose the preferred edge and eliminate 
2-cycles. 

It then repeats the following steps: find a path in the derived graph using 14.71 findpath; 
depending on whether the path found is a simple path or contains a cycle, the algorithm pushes 
flow using H3] pushFlowPath or 14.61 pushFlov^^Cvcle. 

Procedure 14.51 pushFlowPath starts with the first source and moves flow along the path 
without accumulating it at any of the intermediate sources, unless, one of the edges on the path 
is saturated. In this case, the source preceding the edge accumulates some surplus. 

Procedure 14.61 pushFlowCycle is more involved. It uses Procedure 14.51 pushFlowPath 
to transfer the surplus from the initial source to the first source that the algorithm encounters 
in the cyclic part of the path (See figure H]). It then computes the number of pushes of flow 
through the cycles that are required to (i) identify bottleneck edges or (ii) push the surplus to 
zero (See figure ED- Note that an extra last push may be required to transfer the flow from the 
first source in the cycle to the one just before the bottleneck edge as illustrated in figure [6l 

When required, /3j is updated using Procedure 14.41 update and preprocessing is done for 
the next iteration. It also checks if a saturated edge ij which has a small value of "jij prior to 
the increase in /3j has the condition that ^jij = Cij — pijfij — < 0, i.e. the edge is on the verge 
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(i.e. 7ij is close to 0) of becoming unsaturated and allows for the edge to become unsaturated 
by designating the edge as a back-edge. 

Each source maintains the effective profit Cij —Pij/3j for all the unsaturated edges in a heap. 
Whenever an edge is saturated, it is removed from the heap. When ^ij is set to 0, it is re- 
inserted into the heap. Each of these operations requires O(logm) time. These operations are 
performed when Pj changes, for any j. 

The algorithms terminates when all the complementary slackness conditions are approxi- 
mately met. The algorithm can be analyzed for correctness as well as complexity using the 
same framework as section 12.21 

The correctness is shown by proving that both primal and dual solutions are feasible through- 
out the execution and complementary slackness conditions are approximately met at termina- 
tion. 

The complexity is analyzed via a charging argument. Each iteration of the algorithm is 
charged to one of the events: clearing of surplus at a source, saturation of a sink or saturation 
of an edge. All such charges are then accounted for in the final analysis. 

The details of the algorithm and the proofs are presented in the following section. 

4.1 Detailed Algorithm 

The modified algorithm for BTS is presented below in detail, followed by proof of correctness 
and complexity analysis. Since BTP is a special case of BTS, the algorithm is applicable to 
both the problems. 

Algorithm 4.1 Modified Auction 

1: Initialize 
2: Preprocess 

3: while there exists a source i such that Y2j fij ^ ^« and aj > do 
4: P ^ FindPath(i) 

5: if P contains a cycle C which is not a 2-cycle at the end of P then 
6: pushFlowPath (P') where P' is the portion of P before C 
7: pushFlowCycle (C) 
8; else 

9; if P does not end with a 2-cycle then 
10: pushFlovi^Path(P) 
11: else 

12: //P ends at a 2-cycle iiji 

13: pushFlowPath(P') where P' excludes the last edge iiji 

14: Eliminate back edge jiii and let y^j^ ^ /3j^ 

15: end if 

16: end if 

17: Update /3j 

18: end while 



4.2 Feasibility of the Solutions 

Lemma 4.1 The primal and dual solutions maintained by algorithm \4.1\ Modified Auction 

are feasible throughout its execution. 
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Algorithm 4.2 Initialize 

1: fij for all ij 
2: for all j 

3: Jij for all ij 

4: ai i— maxj Cij for all i 



Algorithm 4.3 Preprocess 

1: for all i such that J2j fij < a, > do 

2: Find an unsaturated ij such that Cij — PijPj is maximized 

3: Make ij the preferred edge for source i 

4: ai ^ max(0, Cij - PijPj) 

5: / /Remove all 2-cycles unless it involves the only back edge 
6: if 3i' such that y^/j < j5j then 

7: Vij ^ Pj 
8: end if 
9: end for 



Algorithm 4.4 /3, update 



for all j do 

/ /Update j3j if all assigned sources are shipping at the higher value 
if Mi such that edge ij is unsaturated and yij = fij then 

/3, ^/3,(l + e) 

for all i such that edge ij is saturated do 
if Cij — PijPj — ai <0 then 

/ /There is an edge with 7^ close to zero, i.e. there is another edge with a larger 

effective profit 

Designate ij as a back-edge 



/ /Raise Pj to e miuj ^ if it is and sink j is saturated 



end if 
end if 
end for 
end if 
Preprocess 
end for 



if Z]j fijPij = and Pj = then 

Pij 



/3j ^ e miuj ^ 



19 



Algorithm 4.5 pushFlowPath(P) 



Ensure: P = {k,ji,i2,j2---k,ji} or P = {hji . ■ .ii,ji,ii+i} 
1: //(p is the amount that is transferred. InitiaUy, it is the surplus at the first source in the 
path 

r = I 

if P ends at ji then 

r = / - 1 
end if 

for 1 to r do 

//The transferable amount is subject to the forward and backward edge capacity con- 
straints 

(P ^ min(0, Ui^j^ - fiukJik+dk^J^) 

fidk ^ fikJk + ^ 
Vikjk ^ Pjk 

I /(j) is multiplied by the fraction of two prices to reflect the proportional 
//flow reduction on back edge ik+ijk 



13: 

14: 
15 
16 
17: 
18 
19 



fik+ljk fik+ljk ~ ^ 

end for 

if r = ^ — 1 then 

(f) ^ min(0, Ui^j^ - f^^j^ 

fkii ^ fkji ^ 
end if 
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Algorithm 4.6 pushFlowCycle(C) 



Require: C = {ii, Ji, ^2,i2^ ■■■,ikT jk^ik+i = ^1} is a cycle in the derived graph 

1: //Compute pq 

2: Pq ^ Pikik/Phh Wz=l 

/ / Compute limiting number of revolutions for each edge 
for Z 1, A: do 
if I > 1 then 

Pi 

„C , Til Viziz 
Pi ^IWp,^^^.^ 

else 

Pi ^ Piiji I Viih 

end if 

K ^ argmax^(X;Lo Shp'i^PeY < Uij - fi.ji) 
Ri ^ ai:gmax^{Yll=oShpf{pey < fn+d) 
end for 

Rmin ^ m.mi{Ri,R'i) 
II Rmin can be 00 
Cm ^ argmin,(i?;) 
4m ^ argmin,(i?i) 



'•min 



2 • ) 

mm' rrnn) 



1 1 change the flow on each edge to simulate pushing of multiple revolutions 
for z ^ 1 to /c do 
if Rmin / 00 then 

fizjz ^ hzjz + Ef5j" sij^^ip^y 

else 

/ / Compute the convergence of the infinite series 

fizjz ^ fizjz + T,T=o^np^^(p&y 

end if 

Vizjz ^ Pjz ^ 

f'iz+ijz ^ fiz+i,jz ~ ^j=o" ^iiP'^^{P&y 

end for 

if Rmin is finite then 

//Send one more cycle of flow ; this will saturate some edge in one revolution 
pushFlowPath(ii . . . ikjkh) 

Let P {ii,Ji, ...i^-i, i;} where the edge iiji was saturated in previous step 
pushFlowPath(P) 
end if 
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Algorithm 4.7 findPath(f) 



1: repeat 

2: if Oj > or source i already belongs to P then 

3; return P 

4: end if 

5: Add source i to P 

6: Let j be the sink such that ij is Preferred{i) 

7: if (3j = or sink j already belongs to P then 

8; return P 

9: end if 

10: Add sink j to P 

11: Let i' be a source such that i'j is a back edge 

12: i ^ i' 

13: until false 



o 
o 

f 

p 



^1 




Figure 4: First we bring the surplus to the first source in the cycle, possibly saturating an edge 
en route 

Proof: Consider the initial primal solution fij = for all fij (line[T]in Procedure l4.2I Initialize). 
This is clearly a feasible primal solution. 

Thereafter, the flow is increased or decreased in Procedures 14.51 pushFlowPath and 14.61 
pushFlowCycle. Let us consider each possible change. 

1. Increase/decrease of flow along a forward edge in line [9] of Procedure 14. 51 pushFlowPat h . 

(a) The source constraint is satisfied as the increase in flow on edge ikjk is at most the 
surplus generated at source in the previous iteration. 

(b) Also, the increase is at most the available edge capacity therefore edge capacity 
constraint is satisfied. 

(c) The total increase in the price of incoming fiow on sink j is Pi,.j^. 4> (line \TU\i and 
the total decrease is Pi. , . ^'''^'^ (h- Therefore, there is no net change and the budget 
constraint remains satisfied and tight. 
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Figure 5: Send the flow around the cycle as many times as possible without saturating any edge 

(d) The decrease on edge ikjk is fline [T3l [T9l) where 6 is at most fi.., 7. ^''=+1-''' . 

Therefore the decrease is at most fi^^^j,.- Thus the positivity constraint on edge 
flows is satisfied. 

2. Increase/Decrease in flow in Procedure 14.61 pushFlowCycle. Everytime the flow is 
pushed, the outgoing flow equals the incoming flow on every source and a proportional 
amount out of every sink. Therefore, not only are the source constraints satisfied, but also 
the budget constraints on the sink. The number of times the flow is pushed in the cycle 
is the minimum of the limiting number of revolutions. Therefore the capacity constraints 
are satisfied. 

The last step uses Procedure 14.51 pushFlowPath, which satisfies the constraints as ex- 
plained above. 

Now consider the dual solution and its feasibility. For all z, Oj is initialized to maxj Cij which 
is clearly larger than maxcjj — pijf3j — 'jij as f3j and 7jj are zero to begin with (lines [2] and [3] in 
Procedure 14.21 Initialize) . 

All the instances in the algorithm where these variables change are enumerated below. 
Assuming that the solution is feasible till these changes are made, we show the solution is still 
feasible after the change. 

1. The update of f3j in Procedure 14.41 0j Update at lines HI and [TT] increases it's value. Since 
any increase in /3j implies a decrease in the value Cij — PijPj — ^ij, ai remains larger than 

Cij — Pij/3j — ^ij. 

2. Change in : ^ij is assigned a value equal to max(0, Ojj — Ciji — Pijl^'j). Therefore, for 
every ij, Oij > Cij — PijPj — "fij by construction. 

3. Note that any change in Oi only occurs when there is a change in some it does not 
change independently. When it changes, it is set to maxjj(0,Cij — Pij(3j) (see line H] in 
algorithm 14.31 preprocess). It is therefore, larger than Cij — Pij(3j — ^ij for all ij as 
7ij > 0. 
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Figure 6: Send it one more time possibly saturating an edge and accumulating surplus at the 
source before it 

□ 

4.3 Complementary Slackness 

We prove the following lemma about the satisfaction of complementary slackness conditions. 
Lemma 4.2 The algorithm terminates with the following conditions satisfied 

Vai>0, Oi-^f^j = (19) 

j 

V/3, >0, bj-Y,Pijhj = (20) 

i 

M-fij > 0, Uij - fij = (21) 
V/ij > 0, \cij - Ui -PijPj - lijl < ecij (22) 



Proof: We consider each condition 

1. Source slackness condition (jl9p . These are the terminating condition for algorithm 14.11 
(line[3|). They are, therefore, satisfied at the end. 

2. Sink slackness condition ()20p . Note that /Sj is initialized to zero (Procedure line [2]). 
Thereafter, /3j is changed only when sink j is saturated (Procedure [331 line lll|) . Also, once 
a sink is saturated, it stays saturated. This can be seen by observing that any decrease 
of flow along any back-edge only takes place as a result of an increase in proportional 
amount of flow along some other edge on this sink. 

3. Edge slackness condition (|21|) . These are satisfied initially as 7ij = 0. Thereafter, as long 
as 'jij is greater than zero, the flow does not decrease on this edge as it is never designated 
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a back edge. Only when is set to zero (Procedure 14.41 line[6|), is the edge ahowed to 
become a back edge. 



4. Flow slackness condition (|22|) . Similar to the proof of lemma 12. 7^ we show that the 
following is satisfied. 



Consider unsaturated edges. Whenever a flow is increased along an edge ij it is a preferred 
edge. In Procedure 14.31 line HI whenever an edge is chosen as a preferred edge, Oj is set 
to Cij — PijPj with being zero as it is an unsaturated edge. The slackness condition is 
therefore, satisfied. Subsequently, /3j can rise by a factor of (1 + e) resulting in a rise of 
at most ecij in effective profit, still satisfying ()23p . Once /3j has changed, ij now becomes 
a back edge. Now, f3j can not change any further unless this back edge has zero flow 
(Procedure 14.41 lines 3,4 and 8). If as a consequence of a rise in Oi is set to zero, the 
inequality (|23|) is true and no further flow will be pushed from source i. f3j cannot rise 
until all flow is pushed back to source i. 

Next consider when the edge is saturated. In this case ^ij is set to Cij —pijf3j — ai, thereby 
satisfying ([23]) • However, flow may be pushed back on the saturated edge such that it 
becomes unsaturated and is set to zero. This happens when (3j rises. The change in 
the quantity Cij — pijfij is no more than ecij] thus |7ij| < ecij just prior to being set to 
zero. Thus condition (|23p is satisfied after the change. Condition (j22p immediately follows 
from condition ()23p . 

Hence we conclude the above lemma. □ 

From the above two lemmas we conclude: 

Theorem 4.1 The algorithm determines a primal solution to the BTS problem such that cijfij > 
(1 — e)OPT where OPT is the optimal solution to the BTS problem. 

Proof: The proof is similar to the proof of lemma 12.91 The value of the dual solution is 



= ^(fli - ^ fij)oii + - '^Pijfij)Pj + "^[{uij - fijhij + fijiai +Pijl3j + 7ij)](25) 



When we subtract the value of the primal solution (^ijfiji from the above, the difference 

is 



yfij > 0, ai< Cij - pijPj - 'jij + ecij 



(23) 




(24) 



3 



3 



Cij Oii Pijj3j ^ij) 



(26) 




Thus, the total absolute difference is at most Ai + A2 + A3 + A4 where 
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A3 = lE/5^-(^^--E^'^^/^^-)l 

j i 

A4 = I ^ (Uij - -fijfij)\ 
ij i 

By arguments similar to lemma 12.91 and using lemma 14.21 we get that the difference is at 
most ecijfij. 

□ 



4.4 Complexity 

To prove the complexity we first consider the effect of the pre-processing step during the algo- 
rithm. In this step the algorithm removes every 2-cycle unless there is no other back-edge on 
the sink j. 

The lemma below follows directly: 

Lemma 4.3 If there is a 2-cycle between a source i and a sink j then there is no back-edge 
incident onto j other than ij. 

The algorithm finds path and cycles and pushes flow on edges of the path or cycles. Each 
operation corresponding to a push of flow on an edge either enables a reduction in surplus and 
changes flow to the capacity of the edge. If the flow on a back edge reduces to zero, this leads to 
a change on the dual variables /3j for some sink j. We attempt to charge the various operation 
in finding paths or cycles and the subsequent processing to changes in the dual variables. 

Lemma 4.4 Each operation in the procedures pushFlowPath, pushFlowCycle, findPath 

can be charged to a rise in the value of /3j , for some j such that each rise in the value of (5j is 
charged 0{n'^) operations . 

Proof: We account for the work required for each push of flow along a path or cycle via a charg- 
ing argument. Procedure findPath requires 0{n) steps. Further, during pushFlowPath 
there are at most n pushes of flow along edges of the path before one of these happens 

1. A 2-cycle is encountered. 

If the 2-cycle is between i and j, source i increases the assigned value yij to f5j. Since ji 
was the only back-edge on j (Lemma 14. 3p . there is a change in j3j. We charge the 0{n) 
pushes to this rise in 13 j. 

2. A source with a = is encountered. This case is analyzed together with the next case. 

3. An unsaturated sink is encountered. 

The previous cases can be considered together with this case as they are similar. In 
both these cases, the surplus from a node travels along the edges of the derived graph 
and reaches a source or a sink where no further push of flow is required. The following 
sub-cases arise: 

(a) The surplus at the starting node disappears and no new surplus is created 

In this case, we charge the at most n pushes of flow along the edges of the path to 
the node at which the surplus disappears. We charge such a node only once for the 
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disappearance of the surplus. Note that any subsequent appearance of surplus will 
be charged to the creation of the surplus and is caused only in the case when the flow 
on a back-edge becomes zero, which is the second case explained below. Considering 
all possible source-sink pairs, there are, therefore, at most pushes charged to the 
nodes before either the surplus at each source is removed or a surplus appears. 

(b) There is a back edge ji such that fij becomes zero and a new surplus at i' is created 
where i'j was the preferred edge for i' . 

In this case we charge the push on edges of the path leading to this back edge, i.e. 
a charge of n, to the creation of the surplus at this back edge. By lemma 12.61 this 
back edge can not re-appear unless (3j changes. We allow an additional charge of n 
to pay for the possible disappearance of the flow from i' . 

There are at most n back-edges on sink j. The value /3j changes when the flow on 
each of these edges is reduced to zero. Therefore, the total charge on all the back- 
edges incident out of a sink is at most 2n^ pushes of flow after which there is a change 
in the value This implies we have 2n^ pushes charged to a rise in /3j. 

(c) There is a forward edge ij which becomes saturated to capacity, i.e. fij = Uij. If 
a forward edge is saturated, it stays saturated until it becomes a back edge when 
'jij = 0. When an edge ij becomes a back edge, by lemma [2T6l this edge cannot occur 
again as a back edge until (3j increases. We can, therefore charge this work to the 
corresponding rise in the value of f3j. The rise in could be charged 0{n?) times 
as in the previous case. Note that as a result of this saturating push, surplus could 
be generated on the source just before the edge ij. We put an additional charge of 
0{n) on this edge to pay for the possible disappearance of this surplus in the future. 

4. A cycle is encountered. In pushFlowCycle, all the calculations for determining the 
limiting number of revolutions can be done in 0(n) time, since the cumulative transfer 
ratio can be determined by traversing the cycle once. As noted above, in all the cases that 
arise when flow is pushed around the cycle, either (i) the surplus that has been pushed 
into the cycle reduces to or (ii) the flow on one of the back edges is reduced to zero or 
(iii) the flow on one of the forward edges reaches capacity. The 0(n) amount of work can 
be charged to either the disappearance of surplus at a source, or a reduction of flow to 
zero on a back-edge ij, or the saturation of a forward edge. This implies an increase in 
13 j, subsequently, in the same way as described for simple paths above. 

From the case analysis above, we conclude the total charge on a rise in /3j is 0(n) from case 
1 and O(n^) from case 2,3 and 4, hence proving the lemma. □ 



Lemma 4.5 The algorithm terminates in 0{e "^{ri^ + n log m)m log ?7) time, where 
t^ = max,j(|^)/emini,(|i). 

Proof: From lemmaH31 there are O(n^) charges per rise in the value j3j. Once a phase is over, 
the change in value of /3j causes sources to update the heaps (one at each source). This takes 
no more than 0(n log m) time. 

The update procedure takes 0(1) amortized time by maintaining a set of values of yij, each 
value being either j3j or j3'y Thus, we require 0(n log m) at the end of a phase to update all the 
data structures. The preprocessing requires n steps. 

Each rise causes the quantity /3j to grow by a factor of 1 + e, starting with eminj(^). 
For each sink tj such that Pj > 0, there exists an i such that si ships flow to tj and thus 
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< Oj = Cij —pijPj. This leads to j3j < Cij/pij < maxy — . There are m different f5j. Therefore 

there can be no more than 0(m log^^^ U) total changes in f3j for all j. Combining this with the 
fact that every rise in /3j is charged 0(n^+n log m) amount of work and log^^^ U = 0{e~^ log U), 
we have the result. □ 



5 Concave, Piecewise Linear Profit 

We now describe how to extend the above algorithm to a profit function which is concave and 
piecewise linear. We use the common edge splitting technique in order to reduce the problem to 
the linear profit function. This transformation is very similar to the well known transformation 
of convex cost mincost flows to linear mincost flows, see yy 

Given an instance / of the problem with concave piecewise linear profits, we map it to an 
instance /' of the capacitated but linear profit version. 

Let the profit function be defined as Cijk € Z for the edge ij and interval k each interval 
being of fixed length say total number of intervals being £. 



Figure 7: The profit function for an edge ij. The slope in an interval k is Cijk- Number of 
intervals is C 

We then convert the problem by splitting each edge into C edges, the kth edge being of 
capacity I and profit Cijk- The price of each is the same as the one for edge ij. 

It is easy to see that a solution to instance /' can be easily converted into a solution for 
instance /'. A solution to I', however, can not always be converted into a solution for /. 
Nevertheless, any solution to /' can always be modified into a solution with greater or same 
profit such that the new solution can be transformed into a solution for I. 

If for a solution fijk the following property holds; /Bzi,Z2 such that zi < Z2 and /jj^i < I 
and fijz2 > 0) then we can convert the solution to one for I. 

If this property does not hold then the profit can be increased/kept the same while eliminat- 
ing such pairs. Transfer min(Z — fijzi,fijz2) from the edge ijz2 to ijzi. The profit Cijz^ > Cijz2 
because of concavity and thus we obtain a more profitable solution. Repeating these transfers, 
we get the property that ^zi, Z2 such that zi < Z2 and fijz^ < I and fijz2 > 0. 

This approach can also be used to map an approximate solution. As a consequence, the 
algorithm UTT] Modified Auction can be used to approximate the above problem after the ap- 
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o o 

Cijk Pij 



An edge in the piecewise linear version 




Mapping to linear version, each edge has a capacity 1 

Figure 8: Construction of multiple edges corresponding to one original edge. 

propriate transformation. The transformation, however, does involve an increase in the number 
of edges. We conclude with the following result. 

Theorem 5.1 The Budgeted Transportation problem where the profit function is piecewise- 
linear and concave, can be approximated within a factor o/(l— e) in 0(e~l(ri^+n log m)>Cm log U) 
time where L is the number of intervals in the profit function. 

6 Conclusion 

We have presented an approximation scheme for Budgeted Transportation problem. We further 
generalize it to the capacitated version, which can then used to solve the problem with a non- 
linear profit version. 

The technique used is a kind of primal-dual mechanism based on Auctions. The variables 
are modified in small steps in order to maintain approximate slackness conditions. We use 
augmenting path based mechanism to improve the complexity of the scheme. However, because 
of the nature of the problem, the cycles may arise in addition to simple paths. We have shown 
that the cycles can be handled without too much work. 

An interesting open question is whether it is possible to solve the problem exactly and as a 
consequence, generalized flow using the above scheme. 
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7 Appendix-Relation to Generalized Flow 

Apart from being a natural extension of the Transportation Problem BTP is also related to 
a well known flow problem, generalized flows. As mentioned above, it is a special case of 
Generalized Flow. It is interesting to see, however, that there also exists a reverse relationship. 
We show how to transform the generalized flow problem to Budgeted Transportation. 



Mincost Generalized Flow: We are given a digraph G{V,A), a cost function c : ^4 — > R, a 
capacity function u : A ^ M"*", a multiplier function jj, : A ^ M"*" a source s G V, the supply at 
the source dg, a sink t € F and the demand dt at the sink. The goal is to find a flow function, 
f : A ^ M"*" such that the flow is conserved at the nodes, is multiplied across the arcs and meets 
the supply and demand constraints, while minimizing the total cost of the flow. This problem 
can be expressed as the following linear program. 



minimize / -^^jj 
subject to : 



Y.i'i^k-Y.fik = ^jeV/{s,t} (27) 

i k 

Y^fsk = ds (28) 
k 

i 

hj < Uij VijGA (30) 
fij > Vijeyl (31) 



Mincost Budgeted Transportation: This problem is a more generalized version of the 
conventional transportation problem. We are given a bipartite graph B, consisting of sources 
5 = {sj} and sinks T = {tj}. A supply function a : 5 — )■ M"*", a budget function bj — >■ M"*", a 
cost function c : 5 x T — )■ R and a price function p : S x T ^ M.. The goal is to come up with 
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a flow function f : S x T ^ such that the total flow going out of a source is equal to the 
supply, the total price of the flow coming into a sink is equal to the budget, and the cost of the 
flow is minimized. The problem is stated below as an LP. 



minimize 



■ijfij 



ieS,j£T 

subject to : 



fij = 0,1 i€ S 



(32) 



^Pijhj = hj jeT 



(33) 



fij>0 ieSjeT 



(34) 



Wc first show that an MCGF problem instance transforms to a minimum Budgeted Trans- 
portation problem instance. The optima of the two instances are the same and the solutions 
themselves can be mapped to each other. The transformation is valid even for approximate 
solutions. 

Let / be an instance of the MCGF problem. We transform it to an instance /' of Mincost 
Budgeted Transportation problem as follows 

• For each node i in V, we have a corresponding source Si E S with capacity Oj = Ylj "^ij- 

• For each arc in A, we have a corresponding sink in T", tij, with budget hij — Uij. 

• We have one additional sink ts with budget equal to the supply dg in /. It is connected 
to the source corresponding to source node s with an edge whose price is 1 and cost is 0. 

• Each sink tij has two incoming edges from s j and Sj , and the cost and price are as follows 



• The capacity of the source corresponding to the sink in I is equal to the demand dt 

Consider a pair of flow functions F and F' for the instances / and J', respectively. Let these 
be related to each other according to the relations defined for all ij pairs in the instance /. 





(35) 
(36) 



Further, for the source Sg corresponding to the node s 
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Lemma 7.1 For a given solution F to the original instance I, F' is a feasible solution to I' of 
the same cost. 

Proof: We show that F' is a feasible solution to I' by showing that it meets the constraints [32] 
and[33l 

For each source Sj, except the source Ss, the total outgoing flow is the sum of flows to the 
sinks corresponding to outgoing and incoming edges on node j, which is, 

= \^ f + V f 

i k 

Using the mapping of solutions, this is equal to 

i k 

k i k 

= lijfc = Uj [using equation [27] 
k 

Therefore, constraint 1321 is met. 

For the source Ss, a similar analysis shows that constraint 1321 is met, since fsk = dg and 

For each sink tij the total price of the incoming flow is 

fsi,tijl f^ij ~^ fsjfyj 
— fij ~l~ ^ij fij — ^ij 

Therefore, the constraints [33] are met. 
Also, the total cost of F' 

'^(^■fsi,tij + (^SjUjf'sjMj) ^ fijf^ij'^ij/f^ij = '^ijfij 
tij ij ij 

which is same as the cost of F. □ 



Lemma 7.2 For a given solution F' to the instance I' , F is a feasible solution to I and is of 
the same cost. 

Proof: To show that the conservation of flow constraint is met, consider the total incoming 
flow on a node j, which is 



'^fJ-ijfij — f'sj, Ui 



There are two sets of edges incident on source sj. The ones corresponding to incoming edges 
on node j, {sj,tij) and others that corresponding to outgoing edges {sj,tjk). The sum of the 
flow on these is equal to aj. So the the above is 
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jk 



f's t 

Since the total price of incoming flow on sink tjk is equal to f'g. ^.^ H — which is equal 
to the budget of the sink tjk = bjk-, the above becomes, 

f 

[bjk - — ) 

k f^jk 

= ~ X] ^Jfc + fjk since bjk = Ujk 

k k 

which is equal to the outgoing flow Ylkfjk- For the case of the source and sink, there are no 
incoming or outgoing edges, respectively. But, 

i i 

thus ensuring the sink demand in / and 

k k k 

Thus the outgoing flow at the source, 

fsk = ds 

k 

since as = Yjk'^sk- 

Now we show that the capacity constraints are met. For each edge (i, j) G A, from constraint 
30, we have 

f's t 

fs.Uj + — = bij = Uij 

f 

l^ij - ''^''^^ - 

^ fij ^ Uij 



The cost of the solutions are mapped exactly as in the proof of lemma 

O □ 

Bey combining the lemmas 17.11 and 17.21 we conclude the following theorem regarding the 
equivalence of the two problems. 
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Theorem 7.1 Given an instance of the MCGF problem I we can construct an instance of 
Budgeted Transportation Problem such that the optimum solutions of each have the same cost. 

Also, note that the solutions themselves can be mapped to each other. Thus, by solving 
the Budgeted Transportation problem we not only determine value of the optimum solution, 
we can also construct the solution itself. 



7.1 Mincost to Maxprofit Budgeted Transportation 

In order to minimize the cost, we could maximize the negative of the cost —Cij. The maxi- 
mization, however, does not clear the sources and as such conservation of flow in the mapped 
solution to the original problem would not be achieved. 

We can get over this problem by minimizing instead, M — Cij where M is a very large 
constant. Since every edge has a very large profit, the optimum will clear all the sources. Any 
error 5 in clearance will produce a negative term M5 in the total profit. By choosing an M 
large enough, 5 can be made negligibly small; i.e. smaller than the granularity of values in the 
solution of a Linear Program. We refer to [21] for the discussion and bound on this granularity. 

The cost of the optimum would be M^-a, — Ylij'^ijfij^ where Ylij^ijfij is an optimum 
solution for the mincost instance. We note that an approximate solution to Maxprofit Budgeted 
Transportation, would not immediately lead to an approximation algorithm for Generalized flow. 
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